Microwave anisotropies from the Galactic halo 
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ABSTRACT 

Models in which a large fraction of the Galactic dark matter takes the form of cold gas 
clouds imply that there is thermal microwave emission from the Galactic dark halo. Such models 
can therefore be directly constrained by data on the microwave sky, and in particular the very 
sensitive observations of microwave anisotropies which are now being made. To this end we 
have computed the anisotropy power-spectrum expected for a Galactic dark halo made of cold, 
dense gas clouds, including the effects of clustering with a CDM-like mass spectrum of mini- 
halo substructure. The power-spectrum displays two peaks: one, at I ~ 50, is the Poisson noise 
for the mini-halos, and the second, much larger and at much higher I, is the Poisson noise of 
the individual clouds. Because it appears on small (milli-arcsecond) angular scales, where the 
instrumental sensitivity is inevitably very poor, the latter signal is not directly detectable. By 
contrast, clusters of cold gas clouds may contribute significantly to the observed anisotropies if 
their emission has a grey-body spectrum. In this case the peak fluctuation, at I ~ 50, amounts 
to 4/| sin 6 1 /iK in the Ray leigh- Jeans limit, and is the dominant Galactic foreground between 40 
and 80 GHz. It will be possible to constrain this foreground component using low-latitude data 
from the MAP satellite, providing that its spectrum conforms to a grey-body. If the spectrum 
is "dusty" there will be relatively little power at frequencies below the thermal peak, and in this 
case the predicted anisotropies are shown to be negligible. 

Subject headings: cosmic microwave background — galaxies: halos — dark matter 



1. Introduction 

Studies of the microwave sky play a crit- 
ical role in cosmology, and measurements of 
the degree-scale microwave anisotropies can pro- 
vide tight constraints on the cosmological model 
(e.g. Hu and Dodelson 2002). It is there- 
fore particularly important that we understand 
the power-spectrum arising from foreground mi- 
crowave sources, and much attention has already 
been given to this topic (e.g. Tegmark and Ef- 
stathiou 1996; de Oliveira-Costa and Tegmark 
1999). The various foreground sources which have 
been considered to date include: normal galaxies 



and Active Galactic Nuclei; Galactic synchrotron 
emission; Galactic free-free emission; and Galac- 
tic dust emission. At frequencies in the range 
50 < v < 150 GHz, none of these sources is 
thought to make a significant contribution to the 
power-spectrum for multipoles I < 500 (Tegmark 
and Efstathiou 1996). This conclusion must be re- 
garded as tentative in the case of the extragalactic 
sources, because the difficulty of undertaking large 
area surveys at high radio-frequencies means that 
the source populations are poorly constrained in 
this frequency range — see Toffolatti et al (1999) 
for a discussion. The Galactic emission, on the 
other hand, is known to be very significant at 
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both low frequencies (synchrotron) and high fre- 
quencies (dust); and, in fact, at all frequencies 
at low Galactic latitudes. However, observing at 
high/low frequencies yields template distributions 
for those contaminants, and these templates can 
be cross-correlated with data in the preferred fre- 
quency band to estimate their contamination level. 
This procedure has provided robust confirmation 
that, at high latitudes, the recognised Galactic 
foregrounds do not significantly contaminate the 
COBE DMR data (Kogut et al 1996). 

Although robust, the method we have just de- 
scribed clearly cannot be applied if there is no 
template for the foreground. This is the circum- 
stance we are faced with if the Galactic dark mat- 
ter has a substantial component in the form of cold 
gas clouds, as has been proposed by a number of 
authors (Pfcnnigcr, Combes and Martinet 1994; 
De Paolis ct al 1995; Henriksen and Widrow 1995; 
Gerhard and Silk 1996; Walker and Wardle 1998; 
Sciama 2000a). Such clouds emit primarily in the 
microwave band, and that emission is thermal at 
temperatures of only a few degrees; consequently 
the magnitude of this putative foreground is at 
present only subject to very weak direct observa- 
tional constraints. Indirect constraints - from the 
small amplitude of fluctuations in the Cosmic Mi- 
crowave Background (CMB), and from Big Bang 
Nucleosynthesis - are generally supposed to ex- 
clude any significant amount of baryonic dark mat- 
ter (e.g. Turner and Tyson 1999), but these argu- 
ments are not entirely free of loopholes (Hogan 
1993; Walker and Wardle 1999) and direct con- 
straints are desirable. Given the great sensitivity 
of the microwave data which are now being ac- 
quired, it is therefore prudent to consider what 
Galactic microwave emission is expected for mod- 
els in which the dark matter is composed of cold 
gas. 

To proceed with a calculation we need to spec- 
ify a particular model. The success of the now- 
standard structure formation paradigm, exempli- 
fied by the Cold Dark Matter model (e.g. Pee- 
bles 1992), argues that any acceptable model must 
possess clustering properties that are similar to 
those of CDM, in order to match the data on 
large-scale-structure, and consequently we utilise 
the properties of CDM halos as a guide. It then 
remains to specify the properties of the individ- 
ual clouds. Radio-wave scintillation data - specif- 



ically the Extreme Scattering Events (Fiedler et al 
1987; Fiedler et al 1994) - suggest a vast popula- 
tion of ~ AU-sized clouds, each of planetary mass 
(Walker and Wardle 1998), and we assume that 
these properties are appropriate to the individual 
clouds. 

The present paper is organised as follows: in §2 
we derive the angular spectrum for an unclustered 
cloud distribution; in §3 the influence of CDM-likc 
clustering is considered; §4 combines the results of 
§§2,3 into a prediction for the root-mean-square 
temperature anisotropy spectrum, under the as- 
sumption of grey-body emission, and considers the 
constraints which existing data place upon the in- 
put model; §5 then considers how these results are 
modified in the case where the clouds have "dusty" 
emission spectra. 

2. Power from individual clouds 

As a tool to describe the statistical properties 
of the microwave sky, it is conventional to employ 
the power-spectrum 
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where the ai m are coefficients in the expansion of 
the brightness temperature field, 5T, in terms of 
the spherical harmonic functions Yi m : 

oo I 

5T(M) = £ E aimYi m (0,<f>). (2) 

i=0 m=-l 

We will model only the power-spectrum which 
would be derived from measurements over a small 
patch of sky, so that we need only consider the 
case I ^> 1, for which some mathematical simplifi- 
cations are possible. The restriction to large / also 
means that we do not have to model the struc- 
ture of the Galaxy in any detailed way, because 
we are considering angular scales which are suffi- 
ciently small that the source (cloud) positions are 
essentially random. 

The starting point for our calculation is the 
white-noise spectrum (C; = Cj°, independent of 
I) which is introduced by a random distribution of 
point sources, each of flux F, with a mean number 
per unit solid-angle of dN/dft: 

Cr=F^. (3) 
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If the individual sources are not point-like, but 
each source is axisymmetric with power-spectrum 
F 2 Cf (so that Cf is the power-spectrum for a 
source of unit flux), then the power-spectrum of 
the source population is C/ = 47rCf C™ . Now if 
we consider sources at various distances, D, with 
corresponding number density n, it follows that 
S(dN/dfl) = D 2 nSD, and hence 



^ = 4irCj s D 2 nF 2 . 
dD 1 



(4) 



Integrating along the line-of-sight then yields the 
power-spectrum of this source population. 

The gas clouds we are considering are ex- 
pected to cool primarily through thermal con- 
tinuum emission in the microwave band, with this 
emission arising from precipitates of molecular hy- 
drogen (Wardle and Walker 1999). The bulk of 
the interior of a cloud is too hot for solid hydro- 
gen to exist; all the radiation is expected to come 
from a thin, nearly isothermal atmosphere. There 
is at present no quantitative description of the in- 
ternal cloud structure, and consequently we have 
adopted a simple model of the surface brightness 
profile, namely a uniform intensity disk. Provided 
that the cloud radius is R <C D, this model yields 
a power-spectrum 



C s_l_ 
Cl "47T 



2J 1 (IR /D) 
IRo/D 



(5) 



for a single cloud, with Ji being a Bessel func- 
tion of the first kind. For / ^> D/R Q this model 
is a poor one in that it over-predicts the power 
spectrum — a result of the discontinuity in the 
intensity profile at the limb of the cloud. To coun- 
teract this tendency to overpredict the power at 
large I, we have simply truncated the power spec- 
trum at IR = 10D. This ad hoc treatment of 
the high frequency behaviour means that the pre- 
dicted power-spectrum for the cloud population 
is unreliable in the regime IR 3> (D); however 
the corresponding angular scales are so small (sub- 
milli-arcsecond; see §4) that this limitation is not 
important in practice. 

In equilibrium the luminosity of any cloud must 
be balanced by heating, so a cloud of mass M 
satisfies TM a = 4ttD 2 F, where T is the heating 
rate per unit mass. This heating is expected to be 
due primarily to cosmic-rays (Wardle and Walker 
1999; Sciama 2000a). In this paper we consider 



only lines-of-sight at high or mid Galactic lati- 
tude, and for these directions it is reasonable to 
model the cosmic-ray density variations with a 
simple analytic function; we employ the form T = 
ro exp(— \z\/h), where z = D sin b is the height 
above the Galactic disk at Galactic latitude b. A 
more realistic model is offered by a sum of thin 
and thick disks (Higdon 1979; Bcuermann, Kan- 
bach and Berkhuijsen 1985), but for our illustra- 
tive calculations this level of detail is unwarranted 
and we have adopted a single component model, 
with a scale-height of h = 3 kpc (corresponding 
to the thick disk). This choice of scale-height is a 
compromise between the (larger) values deduced 
from analysis of the radio synchrotron data (Hig- 
don 1979; Beuermann, Kanbach and Berkhuijsen 
1985), and the (smaller) values inferred from exist- 
ing cosmic-ray propagation models (Webber, Lee 
and Gupta 1992). Note, however, that if the puta- 
tive cloud population does indeed exist, then these 
models of cosmic-ray propagation need to be re- 
vised (Sciama 2000a, b). 

To complete our prescription we need to spec- 
ify the number density of the clouds. For this 
purpose we employ the collisional isothermal halo 
model of Walker (1999: W99). In this model, or 
indeed any quasi-spherical halo model, the density 
of dark matter varies only slowly in comparison to 
the scale-height of the cosmic-ray disk, and this 
variation can be neglected: 
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(6) 



where V ~ 220 km s _1 is the rotation speed for the 
Galactic disk, r gc ~ 8.5 kpc is the distance to the 
Galactic Centre, and r<2 is the core radius of the 
dark halo. For the Galaxy, the latter is specified 
as rd — 6.2 kpc in W99's model. 

For a homogeneous source population, with 
n ~ constant, and F oc 1/D 2 , equation 4 in- 
dicates that the low-frequency power (for which 
47rC* ~ 1) is dominated by the nearest clouds 
within the field-of-view, Afi. The lower limit on 
the integration over D must therefore not be taken 
as zero, but as D min , the expected distance of the 
closest cloud, with 
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(7) 



For a field of Ail <~ 0.1 sr, similar to the coverage 
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of the BOOMERanG experiment (de Bernardis et 
al 2000), for example, we have D m i n ~ 0.57 pc. 

Using the descriptions given above, we can com- 
pute the power-spectrum for our model of unclus- 
tered gas clouds, Cf: 
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2J 1 (k/x) 
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where x = 2D/h, so that y = 2D min /h, and k 
2lR a /h. The power level at zero frequency is 



(8) 
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(9) 



Figure 1 shows a plot of SI = k^Cf{k)/Cf{Q), 
the standard-deviation of the intensity fluctua- 
tions. In constructing this figure we have taken 
| sin b | = 1 , and we have adopted the parameters of 
the model given above, for which D min ~ 0.57 pc, 
hence y ~ 3.8 x 10~ 4 . 

Three distinct regimes of behaviour can be seen 
in figure 1. For k -C y, all of the clouds are 
expected to be unresolved on the corresponding 
angular scales, so the power-spectrum is simply 
C"(fc) ~ C"(0) and the r.m.s. intensity fluctua- 
tions are proportional to k. For k <~ y some of 
the closest clouds are resolved out, and for y <C 
k < 1/| sin6| the power flattens to (k) / (0) ~ 
2y/nk, so the r.m.s. intensity is SI cx \fk. It 
should be noted that up to this point the results 
are independent of the scale-height, h, and the 
line-of-sight, i.e. the latitude 6; in other words 
the behaviour is completely determined by the lo- 
cal conditions (dark matter density and cosmic- 
ray heating rate). However, the location (k = 
k p ~ -v/2/l sin 6|) and height of the peak fluctu- 
ations are influenced by h and b: clouds which lie 
outside the cosmic-ray disk have very low luminos- 
ity and consequently contribute little to the tem- 
perature anisotropies, so on angular scales smaller 
than those characteristic of a cloud at distance 
D <~ h/\ sin&| the fluctuations decline rapidly. The 
peak value of the anisotropy spectrum can be es- 
timated from C?(k p )/C?{Q) ~ \/2y\ sin6|/7r. If 
we make numerical estimates of these values we 
find that the peak occurs on milli-arcsecond scales 
(I ~ 3 x 10 8 ) and reaches a value ~ 4 mK, as 
discussed in §4. 



2.1. Source counts 

In connecting the foregoing model with obser- 
vations we must take a little care, in that low- 
frequency power is dominated by the brightest 
sources within the field-of-view. If observations 
are made with a sensitive instrument that has 
good angular resolution, then these will be recog- 
nised as point sources and subtracted from the 
map before the anisotropy power-spectrum is es- 
timated. In order to quantify this possibility we 
need to estimate the source counts. For bright 
sources, with D <C h/\smb\, the source luminos- 
ity and density arc uniform over the population, 
hence the usual Euclidean result applies: 
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(10) 



where F max is the expected flux of the brightest 
source on the sky. For the case we are considering, 
the bolometric luminosity is TM Q , so the bright- 
est source is expected to have a flux of F max = 
TM /A , KD^ nin . This result breaks down at low flux 
levels, where the sources lie outside the cosmic-ray 
disk, but F max corresponds to a source which is 
only about 0.1 pc distant, so with a disk scale- 
height of h = 3 kpc the form given in equation 
10 is a good approximation for F 3> 10~ 9 F max . 
Numerical estimates for F max are given in §4. 

2.2. Mean sky brightness 

Using the foregoing model, it is straightforward 
to calculate the mean sky brightness contributed 
by the cloud population, the result is: 



(I) = 



r U 2 /icsc|6| 



r 2 d y 



(11) 



The mean sky brightness thus does not depend ex- 
plicitly on the properties of the individual clouds, 
although there is an implicit dependence of To on 
the column-density of the clouds, as we shall see in 
§4.1. Consequently the same result for (I) is ob- 
tained for the model described in the next section, 
in which clustering of the clouds is considered. 
Using the cosmic-ray heating rate determined in 
§4.1, and an atmospheric temperature of 4.2 K for 
the clouds (§4.2), we deduce a mean sky bright- 
ness temperature, in the Rayleigh-Jeans regime, 
of (Tb) ~ 30 £iK at the Galactic pole. 
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Fig. 1 . — The root-mean-square intensity fluctuations, SI, due to a randomly distributed population of cold 
clouds, heated by cosmic-rays, as described by equation 8. This graph is appropriate for observations of a 
field of size Af2 = 0.1 sr near the Galactic pole (\b\ = 90°). The axes are scaled such that k = 2lR a /h, and 

SI = ky/C?(k)/C?(0). 
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Fig. 2. — The root- mean-square intensity fluctuations, 81 = K^Cf(n)/Cf(0), due to a population of mini- 
halos with a CDM-like mass spectrum. This plot is appropriate to a survey region of size Af2 = 0.1 sr near 
the Galactic pole. At low frequencies (k <§C 1) essentially all of the mini-halos are unresolved and SI ~ k, 
reaching a peak of SI ~ 1 at k ~ 1. This peak reflects the characteristic angular size of the closest mini-halos, 
which is of order 4°. 
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3. Power due to clustering 

The Cold Dark Matter model provides an ele- 
gant framework within which we can understand 
the formation of structure in the Universe (Blu- 
menthal et al 1984; Davis et al 1985; Peebles 
1992). The key ingredients of this model are (i) 
a scale-invariant power-spectrum of density fluc- 
tuations, and (ii) dark matter which initially has 
low velocities. Gravity does the rest. Dense gas 
clouds are expected to satisfy the second of these 
requirements (i.e. they are a type of cold dark 
matter) , provided that they form in the early Uni- 
verse and survive to the present time (Hogan 1993; 
Walker and Wardle 1999). The first requirement 
is naturally satisfied for inflationary models of the 
early Universe, independent of the nature of the 
dark matter. Consequently structures forming in 
a Universe dominated by dense gas clouds may re- 
semble those found in CDM models. This is as 
true for galactic and sub-galactic mass-scales as it 
is for cluster and super-cluster scales. The sub- 
structure within CDM halos follows a power-law 
mass spectrum with index approximately equal to 
-2 (Moore et al 1999; Klypin et al 1999), and 
we adopt this result as our model for the cluster- 
ing properties of dense gas clouds comprising the 
Galactic dark halo. Consequently we expect the 
halo to be populated with a number density n m h 
of mini-halos, such that 



dM M 2 log e (M 2 /Mi) 



(12) 



extending over a mass range Mi < M < M^. 
It should be noted that this formulation implies 
that all of the dark matter is in the form of mini- 
halos, and this is apparently at odds with the re- 
sults of simulations, for which only ~ 10% of the 
halo mass appears in this form (e.g. Ghigna et al 
1998). However, the simulations which have been 
undertaken to date are capable of resolving only 
the upper end of the mini-halo mass range which 
we consider here (see §4.3), and such simulations 
have little to say about the existence of lower mass 
halos. Moreover, if the only relevant physics is 
cold gravitational collapse, acting on a power-law 
spectrum of density perturbations, then there is no 
reason to expect any characteristic mass-scale in 
the clustering, until one reaches masses which are 
so small that the discreteness of the dark matter 
particles begins to play a role. At the very least 



the circumstance described by equation (12) can 
be considered as a limiting case: the total density 
in mini-halos cannot exceed the total dark matter 
density. 

Because these aggregates are much smaller than 
their parent halos, their internal structure has 
not yet been thoroughly investigated with simu- 
lations. However, we know that their density pro- 
files are consistent with tidally-truncated isother- 
mal spheres (Moore et al 1999), and we adopt that 
description. In addition to tidal truncation, the 
density distribution must have a central core as 
a result of collisions occurring between the con- 
stituent clouds (W99). We therefore employ a 
model density profile 



Mexp(-r 2 /2r t 2 ) 
(27r) 3 / 2 r t (r 2 +r 2 ) 



(13) 



where r t is the tidal radius and r c the core ra- 
dius of the mini-halo, and we have assumed that 
t c *C r t . From this three-dimensional profile we 
can determine the column-density seen by a dis- 
tant observer. The microwave intensity is pro- 
portional to this column, so the angular power- 
spectrum for a source of unit flux is 



Cf 



exp(-2lr c /D) 
47r + 8; 2 r 2 /D 2 



(14) 



The features of this power-spectrum are as fol- 
lows. For lr t < D, Cf ~ and for D/r t < 
I <C D/r c we have Ci oc l/l 2 , cutting off expo- 
nentially for / > D/r c . Most of the contribution 
to the intensity fluctuations therefore comes from 
the frequency range D/r t <C I <C D/r c , and the 
existence of a core in the density profile is critical 
in determining the power-spectrum at high wave- 
numbers. If there were no core in the halo density 
profile, the intensity fluctuations would be inde- 
pendent of I in the range l^§> D/r t . 

We need to calculate the tidal radius and the 
core radius; both depend on the mass of the mini- 
halo under consideration. To estimate the tidal 
radius we note that the scale-height of the cosmic- 
ray disk, h, is expected to be small in compar- 
ison with r gc , so we can make the approxima- 
tion that all of the mini-halos within the ficld- 
of-view which make a significant contribution to 
the power-spectrum are at the same Galactocen- 
tric radius r — r gc . Now the tidal radius of each 
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mini-halo is determined by the tides at perigalac- 
ticon, and because r gc is an upper limit to the 
galactocentric radius at that point - i.e. most 
mini-halos will be on orbits which take them much 
closer to the centre of the Galaxy - and simula- 
tions of halo formation find that mini-halo orbits 
have quite large eccentricities (Ghigna et al 1998), 
we take r gc /2 as an estimate of the Galactocentric 
radius at pcrigalacticon. We then follow the usual 
procedure (Binney and Tremaine 1987), leading to 



an estimated tidal radius r 3 := aM, where 
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gc 
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8.0 x 10 21 cm 3 g" 1 . 



(15) 



The core radius of each halo is determined by 
the collisional nature of the dark matter (W99), 
with the velocity dispersion of each mini-halo be- 
ing related to its mass by comparing equation 
13 (for r <C r t ) with the standard form for a 
self-gravitating isothermal sphere. This procedure 
leads to the estimate r 2 = 7M, where 



7 : 



12V2Vt H 
(2tt) 3 / 4 ttS> ; 



~ 2.5 cm 2 g _1 , 



(16) 



with tn — 10 10 years being the elapsed time 
since the mini-halo virialised, and E = M /-kR 2 is 
the mean column-density of the individual clouds. 
W99's model gives E ~ 140 gem -2 , and we have 
used that value in estimating 7. 

The power-spectrum for the whole population 
of mini-halos can be obtained by integrating 



d 2 d 
dD dM 



2 dn m h 2 
dM ' 



(17) 



with respect to M and D. The integration over 
distance is taken out to infinity but, as for the case 
of the individual clouds, the lower limit cannot be 
taken as zero because at low frequencies the power- 
spectrum is dominated by contributions from the 
nearest sources. We can estimate the lower dis- 
tance limit, as a function of mini-halo mass, by 
analogy with equation 7, making the replacement 
n 1 ^ log e (M 2 /Mi)dn m ^/dlog e M, whence 



AO 



^ i „log e (M 2 /M 1 ) 



dn-r, 



1 



d log e M 2 



(18) 



This equation states that the closest mini-halo 
within the mass range Mi to M 2 has a 50% chance 
of being located within a distance D min . If we 



now define z = 2D m i n /h, where h = h/\smb\ 
then with reference to the expected mini-halo mass 
spectrum (equation 12) we can write z 3 := 0M, 
where 
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h 3 Aflp 



2.7 x 10 



-40 



(19) 



In estimating the numerical value of 0, here, we 
have assumed a field-of-view of AO = 0.1 sr at 
the Galactic pole (| sin b | = 1), with the mass lim- 
its Mi ; 2 as given in §4. An important point can 
now be made, even before we calculate the power- 
spectrum of the population. Because the tidal ra- 
dius, rt, and the distance to the closest mini- halo, 
D m in, both vary as M 1 / 3 , the angular radius of the 
nearest mini-halo is r t /D min = {2/h)(a/ 0) 1 / 3 <~ 
4°, independent of M , (and also independent of h 
and b). Remembering that the power-spectrum is 
dominated by contributions from the closest mini- 
halos, we thus expect a peak in the intensity fluc- 
tuations on degree scales. 

We can now integrate equation 17 over D and 
M to obtain our power-spectrum estimate; the re- 
sult is 

(20) 



C? = C?(0)1(k,k c ), 



where 



cf(o) 



and 
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61og e (M 2 /Mi 



■M, 



dx 



exp(- 



G(r 2 gc +r 2 d ) 



(21) 



r 3/2 



/k c x) 



X 2 + K 2 Z 2 



(22) 

having introduced x = 2D/h. The wavenumber 
is now expressed as n = (l/h)(a/0) 1 / 3 (8/n) 1 / 2 , 
and the parameter n c is n c :— a 1 / 3 /3 1 / 6 / \f2Tvy. 
The limits of integration over z are determined 
by the limits of the mini-halo mass-spectrum, so 
z\ 2 := 0M\ t 2, and we have adopted the values of 
Mi ,2 as estimated in §4, leading to z\ = 3.8 x 10~ 3 
and z 2 = 9.1. 

It is easy to verify that I — > 1 as n — > 0. Unfor- 
tunately this low frequency regime (i.e. n <C 1) is 
the only one for which we have been able to deter- 
mine a simple analytic form for J. It is, however, 
straightforward to evaluate 1 numerically, once k c 
is specified. Using the numerical values for the 
various parameters which we have already given 
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(see also the following section), we find k c ~ f .3, 
leading to the fluctuation spectrum shown in fig- 
ure 2. At low frequencies, where essentially all of 
the mini-halos are unresolved, this spectrum dis- 
plays the usual form for Poisson noise from point- 
like sources: 51 oc k. The spectrum then peaks at 
k ~ 1, corresponding to the characteristic angular 
scale of the nearest mini-halos, and declines slowly 
as k increases. 

4. Model parameters 

The anisotropy models we have just presented 
do not incorporate any free parameters: all of the 
various inputs can be estimated - albeit rather 
crudely in some cases - from other types of data. 
We now turn to the business of estimating those 
quantities which have not yet been specified, 
in order to arrive at numerical results for the 
brightness-temperature anisotropy spectrum. 

4.1. Cosmic-ray heating 

For high column-density clouds, such as we are 
considering here (S ~ 140 gem -2 , W99), there is 
significant attenuation of the cosmic-ray flux from 
surface to centre, and so the cosmic-ray heating 
rate differs from that deduced for diffuse inter- 
stellar gas (Umebayashi and Nakano 1981; Sciama 
2000a). We have estimated the cosmic-ray heat- 
ing rate local to the Sun, T , appropriate to dense 
clouds, using the Monte Carlo simulation pro- 
gram GEANT (Ohishi, Mori and Walker 2003; 
http : //wwwinf o . cern. ch/asd/geant). GEANT 
was developed for studying interactions in labo- 
ratory particle detectors, but can be used for as- 
trophysical calculations without modification. By 
using GEANT we are able to study the transport 
of cosmic rays into high column-density clouds, 
with all relevant particle interactions included. 
The main disadvantage of GEANT, for our pur- 
poses, is that it does not conserve energy very 
accurately in the hadronic interactions (GEANT 
Bug Reports #171 & #389). There is a "patch" 
available for the relevant piece of code, GHEISHA 
which treats hadronic interactions, but this patch 
appears to introduce further problems of its own 
(GEANT Bug Report #415), and we have there- 
fore used GEANT with the un-patched version 
of GHEISHA. We can easily calculate the heat- 
ing rate in the limiting case of very dense clouds, 



without using GEANT, because all of the incident 
power is absorbed in this case. We find that in 
this limit GEANT over-estimates the power depo- 
sition by ~25%, and this estimate can be taken 
as a gauge of the accuracy of our computations of 
the heating rate, T. 

For the purposes of estimating the microwave 
emissivity of the gas, we need only calculate the 
total power deposited by the incident cosmic-rays; 
however, the interaction between the cosmic-rays 
and the gas also gives rise to gamma-rays (prin- 
cipally via electron bremmstrahlung and n° pro- 
duction) , and we will report elsewhere on the con- 
sequences of these processes for the interpretation 
of the mean gamma-ray spectrum of diffuse Galac- 
tic emission (Ohishi, Mori and Walker 2003), and 
for the interpretation of the unidentified EGRET 
sources (Walker, Mori and Ohishi 2003). In our 
calculations of the gas heating rate we approxi- 
mated each cloud as a sphere of uniform density, 
and used GEANT to compute the power deposited 
as a function of the gas density (the radius of each 
cloud being fixed at 1 AU). We adopted the "me- 
dian" cosmic-ray proton spectrum quoted by Mori 
(1997; but note that the units on his equation [3] 
should read cm _2 s -1 sr~ 1 GeV~ 1 ), including ener- 
gies down to 100 MeV, and the electron spectrum 
from Skibo and Ramaty (1993), with a low-energy 
cutoff of 10 MeV; more details of these simula- 
tions will be presented elsewhere (Ohishi, Mori 
and Walker 2003). These cosmic-ray spectra are 
appropriate to the vicinity of the Sun, and it is im- 
plicit in our model that the spectra elsewhere in 
the Galaxy have the same shape, and differ only 
in normalisation. We caution that the cosmic-ray 
spectra which we are employing are not identical 
to the spectra which are directly observed, but are 
"demodulated" , to correct for the influence of the 
solar magnetosphere. At low particle energies the 
demodulation corrections are large, with large as- 
sociated uncertainties. 

Our results are shown in figure 3, in the form 
of heating rate per unit mass, as a function of 
the mean column-density, i.e. To(S). Electron 
and proton contributions are shown separately 
(points) and in combination (solid line). Taken 
separately the electrons and protons display the 
same qualitative behaviour: at low columns the 
specific heating rate is approximately constant, 
while at high columns it declines according to 
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Fig. 3. — The cosmic-ray heating rate for dense gas clouds, as a function of the column-density of the cloud. 
The points are the results of computations with the GEANT Monte Carlo simulation package: diamonds 
show the electron contribution (primary cosmic-ray electrons of energy greater than 10 MeV); and triangles 
show the proton contribution (protons of energy greater than 100 MeV). The sum of these two components 
gives us the total heating rate shown by the crosses (connected with a solid line) . The dashed line shows the 
limiting case, appropriate to very dense gas clouds, in which all of the incident cosmic-ray power is absorbed. 
The fact that the GEANT results lie slightly above this line, for high gas columns, is a reflection of energy 
non-conservation in hadronic interactions within GEANT. 



Tq oc X -1 . For low column-densities our cal- 
culations yield cosmic-ray heating rates of order 
3 x 10~ 4 ergs^ 1 g _1 which are consistent with 
previously reported cosmic-ray ionisation rates ap- 
propriate to diffuse interstellar gas (e.g. Webber 
1998), given an absorbed energy of approximately 
7 eV per ionisation (Glassgold and Langer 1973; 
Cravens and Dalgarno 1978). However, we cau- 
tion that our calculated heating rates in the low 
column-density regime are sensitive to the low- 
energy cutoff of the primary cosmic-ray electrons 
included in the simulation. More generally, this 
finding implies that at low column-density the 
heating rate is sensitive to the form of the adopted 
cosmic-ray electron spectrum at low energies, and 
in this regime the spectrum is poorly constrained. 
The uncertainties in cosmic-ray spectra have been 
highlighted by the recent determination (McCall 
ct al 2003) of an ionisation rate, in a diffuse cloud 
(column ~ 3 x 10~ 3 gem -2 ), which is 40 times 
higher than previously assumed. This suggests 
that T may be as large as 10~ 2 ergs~ 1 g _1 at 
the extreme left of figure 3. However, for the col- 
umn densities which concern us here the heating 
rate is determined primarily by the total energy 
density in cosmic-rays (see below), and it is not 
clear whether the new data imply any significant 



revision of this quantity. 

The behaviour of the heating rate at high/low 
columns is straightforward to understand: for 
small columns the cloud is "optically thin" to the 
cosmic rays, so that the absorbed power is directly 
proportional to the mass of the cloud and the spe- 
cific heating rate is independent of the column. At 
large columns essentially all of the power incident 
on the cloud surface is absorbed, so that the total 
heating rate is proportional to surface area, inde- 
pendent of the cloud mass, and the specific heat- 
ing rate is therefore inversely proportional to the 
column-density. In this regime the computed heat- 
ing rate is not sensitive to the details of the cosmic- 
ray spectrum, being dependent only on the total 
energy density resident in the cosmic rays. We are 
primarily interested in high column-density gas, 
for which the heating is mainly due to protons — 
most of the cosmic-ray energy-density being res- 
ident in protons of energy ~ few GeV. For our 
preferred column-density of £ ~ 140 gem" 2 (see 
§4.2), the heating rate is r ~ 10~ 4 ergs _1 g _1 ; 
80% of this power is contributed by the protons. 
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4.2. Properties of individual clouds 

Individually the cloud mass and radius are 
poorly known, but the mean column-density, S = 
M /ttR 2 , is quite tightly constrained because it 
determines the rate of conversion of dark to vis- 
ible matter (W99). Matching the model to the 
Tully-Fishcr relation for late-type spirals yields 
E ~ 140 gcnT 2 (W99). Essentially all of the 
uncertainty is then in the choice of cloud mass. 
The viable mass range is currently rather large, 
perhaps as broad as 10~ 6 M Q < M a < 1O~ 2 M 
(Wardle and Walker 1999). At present the best 
estimate we can make for the cloud mass comes 
from the optical variability of quasars: Hawkins 
(1993, 1996) pointed out that the observed op- 
tical variations, on a time-scale of years, could 
be interpreted in terms of a near-critical Univer- 
sal density of planetary-mass gravitational lenses. 
In the context of our model, these lenses are 
to be identified with the individual gas clouds 
(whereas Hawkins [1993] suggested that they might 
be primordial black holes). Hawkins's (1993) es- 
timate of the mass of the individual lenses is 
10~ 3 M Q . Schneider (1993), however, has mod- 
elled the same data set and his analysis suggests 
a mass of 10 -4 M if the observed variability is 
interpreted as lensing. Minty (2001) re-examined 
this issue and demonstrated consistency between 
lensing models and the data for lens masses in 
the range 10~ 5 — 10 _4 Mq. The lower end of this 
range is, however, excluded by data on Galactic 
microlensing (Rafikov and Draine 2001), because 
of refraction in the neutral gas (Draine 1998). We 
adopt the value 10~ 4 M Q . This choice then de- 
termines the cloud radius, because the column- 
density is fixed; correspondingly R a ~ 1.4 AU, 
close to our adopted value (§4.1) of 1 AU. 

The other key aspect of the individual clouds 
is their emission spectrum. The model of Wardle 
and Walker (1999), which focusses on the ther- 
mal stability of the clouds, requires that they 
be cooled by thermal continuum emission from 
H 2 precipitates. Here we assume, initially, that 
the hydrogen "snowflakes" are large compared 
to the ~ mm wavelengths of their thermal emis- 
sion, and that the individual snowflakes have a 
large optical depth to absorption at these wave- 
lengths. These assumptions imply grey-body 
emission (/„ oc B v ) at the temperature of the 
"snowflake". (Later, in §5, we will examine how 



the model predictions are modified if we relax 
these assumptions.) The latter is, to a very good 
approximation, just the temperature of the atmo- 
sphere of the cloud. We noted earlier that the 
interior of the cloud is likely to be convective, 
and we can therefore make an estimate of the 
the atmospheric temperature by calculating where 
an n = 3/2 polytropic model crosses the phase- 
equilibrium curve for hydrogen. Fixing the mean 
column density at S = 140 gem -2 , this method 
yields the following estimated atmospheric tem- 
peratures: 5.8, 4.9, 4.2, 3.7, 3.4 K, for clouds 
of mass 10~ 6 , 10~ 5 , 10" 4 , 10~ 3 , 10~ 2 M Q , re- 
spectively. For our preferred cloud mass, then, 
the atmospheric temperature is roughly 4.2 K, 
and this corresponds to a conversion factor of 
ST /SI = 752 Kerg _1 scm 2 sr between brightness- 
temperature, in the Raylcigh- Jeans limit, and 
bolomctric intensity. In the formulation we have 
given, the root-mean-squarc intensity spectrum - 
i.e. the quantity SI = + 1) Ci/2n - has the 
units of bolometric intensity, and must be mul- 
tiplied by this conversion factor to deduce the 
temperature anisotropy. 

Similarly, when considering source counts, the 
conversion factor between flux density and bolo- 
metric flux for our model parameters is F v /F — 
2.3 x 10~ 16 ^ 2 Hz -1 , in the Raylcigh- Jeans regime, 
with v in GHz. Using these conversion factors 
we can estimate the flux of the closest single 
cloud, it is F max ~ 10 nKsr, implying that it 
may just be detectable by the MAP satellite 
(http://map.gsfc.nasa.gov). Sciama (2000a) 
has previously noted that the proposed Planck 
satellite (http : //astro . estec . esa . nl/Planck) 
would easily detect the closest examples of the 
posited clouds. In more conventional units the ex- 
pected flux of the brightest source can be writen 
as F max ~ 120 mJy at 20 GHz, and scaling as v 2 
in the Raylcigh- Jeans regime [y -C 150 GHz). 

4.3. Properties of the minihalos 

The mass-spectrum of the mini-halos is given 
in equation 12, but we have not yet specified the 
range [Mi, M 2 ] over which it applies. It is implicit 
in equation 12 that the (dark) density p is entirely 
composed of mini-halos with masses in the range 
Mi < M < M 2 . The appropriate value of Mi 
is difficult to estimate. To form a stable cluster 
there should be a large number of clouds in each 
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Fig. 4. — Microwave anisotropy spectra due to emission from cold, Galactic clouds. This plot shows the root- 
mean-square brightness-temperature fluctuations, ST = [1(1 + 1)(C" + Cf)/2ir] 1 / 2 , as a function of spherical 
harmonic index, I, in the Raylcigh-Jcans limit. The large peak at I <~ 3 x 10 s is the Poisson noise from 
individual clouds, while the smaller peak at lower I is due to bound clusters of clouds (mini-halos) with an 
assumed CDM-like mass spectrum. The figure shows the predicted spectrum appropriate to a small area of 
sky (Af2 ~ 0.1 sr) at high/low Galactic latitude (solid/dashed lines: |b| = 90°/|6| = 30° respectively), and 
data taken at high latitude with the DASI instrument (diamonds; Pryke et al 2002), and the CBI (triangles; 
Mason et al 2002). Also shown on this panel are the results from Archeops (crosses; Benoit et al 2002), 
based on data covering an area A£l ~ 1 sr above b = 30°. 



mini-halo, although the choice of what constitutes 
a large number is somewhat arbitrary; we have 
chosen M x = 0.1 M Q ~ 10 3 M o . The upper limit is 
easier to estimate: M 2 should be chosen such that 
there is a 50% chance of finding a more massive 
halo within the field-of-view. If the total mass of 
the dark halo of our Galaxy is M to t , which we take 
to be 2 x 10 12 M (Zaritsky 1999), then it follows 



2 4tt M 2 \og e (M 2 /M i y v ; 

and solving yields M 2 ~ 1.4 x 10 9 M Q for A£l = 
0.1 sr. Fortunately our final results are not 
very sensitive to the particular values of Mi j2 
adopted, as they enter principally through the fac- 
tor log e (M 2 /Mi) ~ 23, and secondarily through 
the limits of integration (21,2 oc M^ 2 3 ). 

4.4. Temperature anisotropics and con- 
straints 

With the above numerical estimates, and the 
power-spectrum formulations derived in §§2,3, we 
are now in a position to quantify the temperature 
anisotropics expected in the present model. The 
results are graphed in figure 4, showing both the 



mini-halo contribution, which leads to the peak at 
I ~ 50, and the much larger peak at I ~ 3 x 10 8 is 
the Poisson noise from individual clouds. 

Recalling, from §2.2, that the mean sky bright- 
ness is approximately 30 /iK, in our model, we see 
that the mini-halos introduce fluctuations which 
are small in comparison with the mean intensity, 
while the reverse is true for the fluctuations due 
to the individual clouds. This difference simply re- 
flects the fact that the mini-halos are sufficiently 
large that they cover the entire sky several times 
over, whereas the individual clouds cover only a 
tiny fraction of the sky and the root-mean-square 
intensity is consequently much greater than the 
mean. It is worth emphasising that each of these 
contributions is computed under the assumption 
that all of the dark matter is in the corresponding 
form - i.e. clustered into mini-halos, or entirely 
unclustered - and the calculated contributions are 
in this sense mutually exclusive. However, recog- 
nising that (i) the minihalos are made up of in- 
dividual clouds, (ii) the minihalos are unbiased 
tracers of the mean density, and (iii) the miniha- 
los introduce only a small modulation around the 
mean sky brightness, we see that in practice the 
predicted properties of the ^milli-arcsecond fluc- 
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tuation peak are largely independent of whether or 
not the clouds are clustered into minihalos. If all 
the clouds are clustered into minihalos, then there 
would, however, be many fewer clouds very local 
to the Sun, and so the low-frequency (i.e. small 
I) power from individual clouds would be much 
reduced. 

How do our predictions compare with exist- 
ing data? We do not consider, here, the low- 
order multipoles (I < 10) which were measured by 
COBE (Smoot et al 1992), because the calcula- 
tions we have undertaken are valid only for I ^> 1. 
The COBE detections at small I did, however, 
create great interest in measuring the higher or- 
der multipoles, and after much effort expended on 
special-purpose experiments there are now several 
clear detections of anisotropics on degree scales 
(e.g. de Bernardis et al 2000; Hanany et al 2000; 
Pryke et al 2002; Bcnoit et al 2002). The peak of 
the observed signal - 5T b ~ 70 /xK at I ~ 200 - 
is roughly 20 times larger than the model predic- 
tions at the same angular scale, for high-latitude 
fields, and for many purposes this small contribu- 
tion to the observed power (about 0.3%) can be 
neglected. However, it may still be possible to de- 
tect the predicted mini-halo anisotropics against 
the background of the dominant CMB fluctua- 
tions, if we select the data appropriately. In par- 
ticular the CMB anisotropics fall to ~ 30 /jK at 
multipoles I <~ 50, near the peak of the predicted 
spectrum, while the latter attains a value ~ 8 /xK 
at |6| = 30°. Data collected by the MAP satellite 1 
should be able to reveal this Galactic foreground, 
provided that the spectrum is close to the assumed 
grey-body form (see §5 and figure 5). 

To date there has been only limited interest 
in the arcminute-scale anisotropies, because Silk 
damping is expected to strongly supress any pri- 
mary cosmological anisotropy on these scales (Silk 
1968). Consequently there is at present only a lim- 
ited amount of data in the region 10 3 < I < 10 4 
(Subrahmanyan et al 1997; Dawson et al 2001; 
Mason et al 2002). Here again the observed signal 
is large in comparison with the model: 25 fiK at 
I ~ 2500 (Mason et al 2002); roughly an order of 
magnitude larger than the anisotropies predicted 
by the present model. Consequently data at these 
very high multipoles do not yet provide strong con- 

http: //map. gsfc .nasa.gov 



straints on the model we have presented. More- 
over these angular scales are below the resolution 
of the instrumentation carried by MAP, so there 
is no immediate prospect of a major improvement 
in sensitivity. It would, however, be useful to ob- 
tain further measurements of the power-spectrum 
using ground-based interferometers, particularly 
with a view to constraining its latitude depen- 
dence, because the origin of the observed high- 
frequency power (Mason et al 2002) is not clear at 
present and it remains possible that it is a Galactic 
signal. 

For I » 10 4 the anisotropies predicted by our 
model are dominated by the individual clouds, and 
more specifically the closest examples within the 
area under study; consequently for these angular 
scales the model predicts very little latitude de- 
pendence in the power-spectrum. The domain of 
large / is also technically challenging, in that the 
sensitivity of an interferometer to surface bright- 
ness fluctuations worsens in proportion to I. Thus 
the predicted milli-Kclvin peak on milli-arcsecond 
scales is not immediately open to experimental 
scrutiny. Fortunately the model predictions can 
be put to the test in another way: rather than 
attempting to measure the anisotropies per se we 
can simply study the bright source counts. At low 
frequencies the flux from each cloud increases as 
v 2 , or faster, whereas the flux from non-thermal 
(synchrotron) radio sources typically declines with 
increasing frequency. Thus in a sky-survey for 
compact sources at high frequencies the predicted 
population should stand out. For the purposes of 
testing the model, the ideal approach is to make a 
survey of a large fraction of the sky, with sufficient 
sensitivity to detect sources having flux much less 
than F max . The model then predicts that many 
sources with thermal spectra should be detected, 
and any such sources can be selected out and sub- 
jected to detailed studies to determine their na- 
ture. 

To date surveys at high radio frequencies have 
covered only a very small fraction of the sky with 
the necessary sensitivity — this is simply a con- 
sequence of the rapid decrease in area of the pri- 
mary beam of an antenna, as the observing fre- 
quency increases, requiring a much larger number 
of pointings to cover a given area of sky. For ex- 
ample, the survey of Taylor et al (2001), made 
with the Ryle Telescope, covered only 63 deg 2 at 
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Fig. 5. — Frequency dependence of the root-mean-square temperature anisotropics at Galactic latitude 
|6| = 90° and multipole I ~ 50, for various processes. The curves shown here are for: the microwave 
background anisotropies (dotted line; T = 2.73 K, and 8Tb = 30 /uK); mini-halo anisotropics for a grey-body 
spectrum at T = 4.2 K (solid line); mini-halo anisotropies for a "dusty" spectrum (i.e. /„ oc v A at low 
frequencies) at T = 4.9 K (short-dashed line); emission from dust (T = 21 K) in the diffuse interstellar 
medium (long-dashed line); synchrotron emission (dot-dashed line); and free- free emission (dot-dot-dot- 
dashed line). The last three of these are taken from Tegmark and Efstathiou (1996). 



v = 15 GHz, to a flux limit of 20 mJy Using the 
predicted source-counts for our model (§§2.1,4.2), 
it is straightforward to estimate that the expected 
number of cloud detections for this survey is ap- 
proximately 0.01 and these data therefore have lit- 
tle to say about the present model. A more ex- 
tensive interferometric survey, covering the entire 
southern sky, is planned for later this year. This 
survey is expected to have a flux limit of roughly 
40 mJy at 20 GHz - thus improving substantially 
on the pilot survey (Ekers et al 2003) - and could 
detect examples of the closest clouds if their spec- 
tra are indeed grey bodies (see §5). 

Deep surveys at much higher frequencies have 
also been undertaken, again covering very small 
areas of the sky. Interestingly, Lawrence (2001) 
has pointed out that the "blank field SCUBA 
sources" are consistent with being a very local 
population of Galactic sources, rather than high- 
redshift star-forming galaxies as is usually as- 
sumed. Indeed Lawrence (2001) points out that 
their properties are remarkably similar to those 
expected for the population of cold gas clouds pos- 
tulated by Walker and Wardle (1998); the main 
point of discrepancy noted by Lawrence is that 
he infers a thin disk population, rather than a 
quasi-spherical halo. However, as noted earlier, 
the luminosity of each cloud is determined by the 



cosmic-ray heating rate, and models of the Galac- 
tic synchrotron emission (Higdon 1979; Beuer- 
mann Kanbach and Bcrkhuijsen 1985) suggest 
that the cosmic-ray energy density initially de- 
clines very rapidly (scale-height ~ 0.3 kpc) as one 
moves away from the plane of the disk, changing to 
a more gradual decline further out. Thus the dis- 
crepancy suggested by Lawrence might not be all 
that serious. Clearly problematic, though, in the 
context of the present model, is the high bolomet- 
ric flux of the SCUBA sources. For a 2 mJy source 
density of 10 3 deg -2 - as per Lawrence (2000) - 
assuming Euclidean source counts, the brightest 
source on the sky should be of order 240 Jy. These 
are continuum fluxes at 850 ^m (3.5 x 10 11 Hz), so 
the bolometric flux of the brightest source should 
be of order 10 -9 ergcm _2 s _1 . The present model 
predicts a maximum bolometric flux (§2.1) which 
is only 1% of this value, based on cosmic-ray heat- 
ing. This conclusion is largely independent of 
the spectrum of the emission, within the plausi- 
ble range of possibilities. However, it is important 
to note that Lawrence's analysis excluded grey- 
body spectra for the SCUBA sources, preferring a 
"dusty" thermal emission spectrum (see next sec- 
tion). Lawrence's interpretation of the SCUBA 
data can be tested by surveying a large area of 
the sky, in order to measure the bright source 
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counts, and following up with detailed studies 
of the brightest (presumably the closest) sources. 
The Planck satellite 2 should be able to detect the 
closest of the individual dense gas clouds (Sciama 
2000a). 

5. Discussion 

The calculations presented in §4.4 assume that 
the cold gas emits as a grey-body. This is plau- 
sible, but not certain. Lawrence (2001) has sug- 
gested that the Blank Field SCUBA sources are 
cold, dense gas clouds, with properties similar to 
those which we have considered in this paper, and 
he uses data on these sources to exclude grey-body 
spectra, preferring "dusty" spectra which exhibit 
much less flux at frequencies below the thermal 
peak. In addition to this uncertainty, the atmo- 
spheric temperature could be slightly higher than 
the 4.2 K we have adopted, again diminishing the 
predicted power at low frequencies. We have quan- 
tified these uncertainties by computing the mini- 
halo anisotropics for two models representing the 
plausible range of emitted spectra: a 4.2 K grey- 
body spectrum (the reference model used in the 
foregoing sections of the present paper), and a 
4.9 K "dusty" spectrum (i.e. I u oc v A at low fre- 
quencies) . The results are shown in figure 5 for the 
I ~ 50 peak of the model anisotropies, together 
with models (adapted from Tegmark and Efs- 
tathiou 1996) for the known Galactic foregrounds 
(diffusely distributed dust, synchrotron and free- 
free emission). The plot is made for |b| = 90°, and 
at lower latitudes all of the Galactic contributions 
scale up as 1/| sinfe|. 

From this figure we can see that, in the case of 
a grey-body spectrum, mini-halo anisotropies are 
predicted to be the dominant Galactic foreground 
for observing frequencies in the range 40-80 GHz, 
at multipoles I ~ 50. Furthermore, the predicted 
power contribution amounts to ~ 5% of the cos- 
mic background anisotropies, at \b\ = 30° and 
may therefore be detectable with data from the 
MAP satellite. On the other hand, if the emission 
spectrum of the clouds is "dusty" we can see that 
the predicted mini-halo anisotropies are negligibly 
small at all frequencies, and could not expect to 
be detected by any current or planned instrument. 

2 http : //sci . esa. int/home/planck 



The foregoing considerations are specific to the 
anisotropy power spectrum, but it may be possi- 
ble to detect microwave emission from the mini- 
halos in other ways. In particular we note that 
the brightest mini-halos may be detectable as in- 
dividual sources. In this case those sources should 
also be bright in gamma-rays, and a sensible strat- 
egy is therefore to search for microwave counter- 
parts to the Unidentified EGRET sources (Walker, 
Mori and Ohishi 2003). We also note the pos- 
sibility of detecting the predicted microwave fore- 
ground component via its mean intensity which, as 
noted earlier, is much larger than the anisotropy 
level. For grey-body emission the mean inten- 
sity is roughly 30/ 1 sin b\ fiK in the Rayleigh-Jeans 
regime. 

The fact that the present model naturally yields 
a degree-scale anisotropy peak, in the form of low- 
temperaturc thermal radiation, suggests that it 
might be possible to construct a model in which 
the observed microwave anisotropies are inter- 
preted entirely in terms of emission from mini- 
halos. In support of this idea we note that (i) the 
frequency dependence of the anisotropies would 
be quite small (assuming a grey-body spectrum at 
4.2 K - see figure 5), and (ii) the heating rate of 
the clouds might be rather larger than our esti- 
mate, either because of a very large population of 
low-energy cosmic-rays (see the discussion in §4.1), 
or because there are other heat sources which are 
more important than cosmic-rays. However, a 
generic prediction of any such model is a strong 
latitude dependence of the peak power, and this is 
at odds with existing constraints on any latitude 
dependence (Griffiths and Lineweaver 2003). 

6. Conclusions 

We have calculated the microwave anisotropies 
which arc expected in a model where the dark halo 
of the Galaxy is composed entirely of cold, dense 
gas clouds, heated by cosmic-rays. There are con- 
tributions to the anisotropy arising from the Pois- 
son noise of individual clouds, and from the Pois- 
son noise of clusters of clouds. The latter contri- 
bution has been computed under the assumption 
of a CDM-like mass-spectrum for the clustering, 
and we find that it peaks at I <~ 50. The Pois- 
son noise contributed by individual clouds is, for 
(«3x 10 8 , isotropic and sufficiently small that 
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it will be extremely difficult to detect. A better 
approach to constraining this component is to at- 
tempt to detect the individual clouds themselves, 
and this is best done with high-resolution all-sky 
surveys at high frequencies. The Planck satellite 
will collect data of this type and should detect 
> 10 such sources. 

The predicted power due to clusters of clouds 
increases toward the Galactic plane as 1/ sin 2 b. If 
the clouds emit as grey-bodies, then such clusters 
are predicted to be the dominant Galactic fore- 
ground over the range 40-80 GHz, at multipoles 
I <~ 50. For these circumstances the predicted 
power contribution is ^5% at |6| = 30°, and by 
virtue of its latitude dependence this component 
may be detectable in MAP data. On the other 
hand, if the clouds have a "dusty" emission spec- 
trum, then their contribution to the anisotropies is 
too small to be detected, being swamped by other 
Galactic foregrounds at all frequencies. 

We thank Mark Wardle for providing his ra- 
dial temperature profile for a polytrope with index 
n = 3/2, and Lister Staveley-Smith for a critical 
reading of the manuscript. 
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